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We analyse the effects of intrinsic decoherence on the probability of generating a tripartite GHZ 
state using a cool trapped ion coupled to a single mode of the cavity field and interacting with a 
resonant laser field. Milburn equation is solved for this tripartite system to obtain the time evolution 
of density matrix as a function of cavity-ion coupling, laser-ion coupling, and the size of unitary time 
step relative to the time scale determined by system parameters. Starting with the system prepared 
initially in a separable state, density matrix is used to calculate the probability of tripartite GHZ 
state generation using coupling strengths reported in a recent experiment. 



Interest in multipartite maximally entangled states is motivated by possible use in quantum communication 
quantum dense coding Q and quantum teleportation 0. Advances in experimental techniques of trapping, cooling, 
preparation and manipulation of ionic quantum statesjj, |5|, |6(, through interaction with external lasers, have made 
the generation of maximally entangled states a possibility. Coupling of trapped ion to quantized field inside an optical 
cavity has also been achieved 0. With the ion trap placed inside a high finesse cavity, we have at hand a tripartite 
system with additional control mechanism offered by ion cavity coupling. This is the basic unit of a more complex 
multiple function system in which quantum states of several trapped ions placed coupled to an optical cavity may be 
manipulated in a controlled manner. In a letterQ, we proposed a scheme to generate three qubit maximally entangled 
GHZ state, using trapped ion interacting with a resonant external laser and sideband tuned single mode of a cavity 
field. Experimental implementation of such a proposal is constrained by the effects of quantum decoherence sources. 
Suppression of quantum coherence occurs due to random and unknown perturbations of system Hamiltonian. In a 
GHZ state generation experiment, the average effect of decoherence sources on probability of finding the tripartite 
system in GHZ state is of interest. In this paper, we investigate the effects of intrinsic decoherence on the composite 
system dynamics using Milburn model 9] . Milburn model of intrinsic decoherence is based on the assumption that a 
closed quantum system evolves in a stochastic sequence of identical unitary transformations. The size of the unitary 
step along with the energy spectrum of the quantum system, determine the decoherence rate. The model applied 
by Milburn to free evolution of a closed quantum system, has since been used to study decoherence effects in JCM 
model with a single quantized field mode [ToL ITlj. decoherence in multiphoton JCM model^^j decoherence of open 
systems |T3, decoherence of a two-level particle coupled to photon field^l an d decoherence of quantum gate 0. We 
solve the Milburn equation exactly for a cool trapped ion coupled to a single mode of the cavity field and interacting 
with a resonant laser field. The evolution dynamics of tripartite system is a function of laser-ion coupling, cavity ion 
coupling, and size of unitary step relative to the time scale determined by system parameters. The system is prepared 
initially in a separable state having ionic center of mass and cavity mode in respective number states, while the ion 
is in ground state. From the density matrix containing intrinsic decoherence effects, we calculate the probability of 
tripartite GHZ state generation using coupling strengths reported in recent experiments^- 

Cold ions in a linear trap 0, 0, @ offer a promising physical system for quantum information processing, as each 
ion allows two qubit state manipulation. With an ion trap placed inside a high finesse optical cavity, we obtain 
a three qubit system with additional control mechanism offered by quantized cavity field. The interaction picture 
Hamiltonian for this system is given in a detailed form in Eq. (4) of ref. fl6l ] . On the time scale involving only a few 
cycles of unitary time evolution, the effective dimensions of three subsystems namely ionic internal state (A), state 
of center of mass motion (B), and cavity field state (C) depend not only on the ion type but also on the relevant 
part of the ion-field interaction. Presently, we consider the resonant interaction of a two level cold trapped ion having 
transition frequency too, with an external laser of frequency lul — ujq, while the cavity field is tuned to red sideband 
of ionic vibrational motion, luq — uj c = v. Here v is the frequency of one dimensional harmonic oscillator trapping 
potential. The Lamb-Dicke (LD) parameters relative to the laser field and the cavity field are denoted by tjl and r\ c 
respectively. Working in the Lamb-Dicke regime, t]l <C 1 and rj c <C 1, the dynamics of the trapped two-level ion is 



* shelly@uel.b jshellyi9uel. br| 
tnsharma@uel.bi Your e-mail address 



2 



governed by interaction Hamiltonian, 

Hj = hfl[a + + er_] + hgr] c a + ba + a_w a) 



(1) 



where a' (a) and are creation (annihilation) operators for vibrational phonon and cavity field photon respectively. 
The ion-phonon and ion-cavity coupling constants are ft and g, whereas Uk(k = z, +, — ) arc the Pauli operators 
qualifying the internal state of the ion. 

To obtain time evolution of the system, we work in the basis \g, m, n) , \e,m,n) , \g, m — 1, n — 1) , and 
|e, m— l,n— 1) . Here m,n — 0, l,..,oo denote the state of ionic vibrational motion and quantized cavity field, 
respectively. The matrix representation of interaction operator Hi in the chosen basis, to the lowest order in t\l and 
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The unitary transformation that diagonalizes Hi is easily obtained and yields the eigenvectors, Hi$ p = E p $ p , 
(p = 1,4). Milburn 9] postulates that on sufficiently short time steps, the quantum system evolves stochastically and 
the change in density operator in a time interval (t, t + r) is given by 



p(t + t) = 7 cxp ( - 1 -9(t)Hi\ p{t) exp ( ^{t)Hi ) , 



(3) 



9(t) being some function of r. In standard quantum mechanics 0(t) — r. In Milburn model lim 0(r) — I/7, a constant 

minimum time step with mean frequency 7. The evolution of density operator p(t) in Milburn model , is governed by 
the equation 



<'x,> j -J^hA p(i)exp (^-Hi ) /,(/) 



and the rate of change of density operator to first order in 7 1 as obtained from Eq. is 
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Transforming to the basis in which Hi is diagonal and integrating Eq. JSJ), we get 

p,q=1.4 



where 



(pW) M = {p(0)L o exp 
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The presence of factors of the type (E p — E q /Ti) in the exponential indicates that energy spectrum of interaction 
operator Hi plays an important role in determining the intrinsic decoherence rate. By solving Eq. JfjJ, for the initial 
state p(0) = \g, m — l,n — 1) (g, m — 1, n — 1|, we get 
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TABLE I: Panzf) peak values at T = (tt/4, 3tt/4), for an/7 = 0.001, 0.005, 0.01, 0.1. 
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Pghz 


Pghz 


(= 011/7) 


(in nano sec) 


T = Tl/4: 


T = 3tt/4 


No intrinsic decoherence 1.0 


1.0 


0.001 


0.43 


0.99 


0.94 


0.005 


2.15 


0.94 


0.78 


0.01 


4.32 


0.89 


0.65 


0.1 


43.20 


0.53 


0.37 
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Next we consider a special case with the ion in its ground state occupying the lowest energy trap state, while the 
cavity is prepared in vacuum state. The density operator for the initial state is p(0) = \g, 0, 0) (g, 0, 0| ,(m = 1, n = 1). 
Without taking into consideration the decoherence effects, the interaction time t needed to generate the maximally 
entangled tripartite GHZ state depends on the the ratio a = (/in/an) , involving parameters characterizing the ion- 
laser and ion-cavity interaction strengths. For interaction time t p such that ant p = 5, /tnfp = ptr,p — 1, 2, (a = 4) 
the system is found to be in maximally entangled tripartite two mode GHZ state 
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Using Eq. Q, we can easily write down the density operator p(t) for p(0) = |g,0,0) (p,0,0|. Next we evaluate the 
probability PGHzif) — tr(j>(t)paHz) of finding the system in the state pghz — Y&ghz) (^ghz\ and obtain 
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In the limit an/7 ~~ * corresponding to continuous time evolution, Eq. J7J) and Eq. 10 reduce to decoherence 
free time evolution of density operator and Pghz respectively. Fig. 1 displays Pghz if) as a function of scaled time 
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FIG. 1: PoHz(t) as a function of scaled time variable T(= ant) in degrees, for the choice jUii\<tn = 4, and intrinsic decoherence 
parameter R = an/7 = 0.001,0.005,0.01,0.1. 



variable T(= ant), for the choice /iu\an = 4, and intrinsic decoherence parameter R = an/7 = 0.001, 0.005, 0.01, 0.1. 
Evidently the evolution dynamics of PGHz{t) is sensitive to changes in intrinsic decoherence parameter R. With 
increase in the value of R the peaks in PcHz{t) plots are seen to get lower. Using the value VL = 8.95MHz and 
recalling that /x m „ = y a 1 nn + il 2 , we calculate the mean frequency 7 of minimum time step. Table 1. lists the 
PGHz{i) peak values at T = 7r/4, 37r/4, for the choices an/7 = 0.001,0.005,0.01,0.1. As expected, the presence of 
decay factors in Eq. © causes the PGHz(t) peak values to decrease with increase in I/7. For an/7 = 0.1, decoherence 
effects dominate the scene resulting in a low PGHz(t) peak value of 0.53. We may recall here that for the coupling 
value determined by f2 = 8.95MHz, the angle T ~ ir/A corresponds to an interaction time of t\ ~ 0.34/is. The 
probability of finding the system in GHZ state can be measured by cavity-photon measurement combined with atomic 
population inversion measurement. 



I. CONCLUSIONS 



We have analysed the evolution dynamics of a single two-level trapped ion coupled to an optical cavity interacting 
with a laser field and a single cavity mode. Assuming that the quantum system evolves in a stochastic sequence of 
identical unitary transformations |9J], the probability of generating a maximally entangled tripartite GHZ state has 
been calculated. We have used infraction parameters reported in a recent experiment 0. As expected, the presence of 
decay factors in time evolution equation, causes the Pqhz (t) peak values to decrease with increase in the average size 
of the minimum time step. The evolution dynamics of PcHz(t) is also sensitive to changes in intrinsic decoherence 
parameter R. Besides that the energy spectrum of system hamiltonian determines the decoherence rate for a given 
initial state preparation. For the case at hand, that is ion initially in its ground state occupying the lowest energy 
trap level, while the cavity is prepared in vacuum state, the ratio R = {\gTj c )/"i is crucial in determining the rate 
at which GHZ state generation probability PGHz(t), decays with time. Even for a large value of aii/7 = 0.01, we 
obtain GHZ state generation probability of 0.89. Decoherence effects due to cavity decay, heating and other random 
sources of decoherence on Pghz have not been considered. 
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